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Luigi Stella^ '^ and Mario Vietri^ 
ABSTRACT 

We show in this Letter that relativistic dragging of inertial frames around fast rotating coUapsed 
stars is substantial and can give rise to observable effects. We apply this to the kHz quasi periodic 
oscillations (QPOs) sources, low mass X-ray binaries (LMXRBs) containing an accreting neutron 
star. Within the beat frequency model, both the Keplerian frequency of the innermost region of the 
' accretion disk (~ 0.3 — 1.2 kHz) and the neutron star spin frequency 0.3 — 0.4 kHz) are directly 

0^ ■ observed. From these the Lense-Thirring precession frequency (tens of Hz) of the same material 

in the innermost disk regions which gives rise to the kHz QPOs is determined within a factor of 
Oh' ~ 4, depending on the neutron star equation of state. The classical contribution from neutron star 

^ oblateness decreases the precession frequency slightly. The broad peaks at frequencies ~ 20 — 40 Hz 

in the power spectra of the "Atoll" -sources 4U 1728-34, 4U 0614-^091 and KS 1731-260 and their 
variations with the higher kHz QPO frequency are well matched by Lense-Thirring precession of 
material in the innermost disk region. We also suggest that the '~ 15 — 50 Hz horizontal branch 
QPOs of GX 5-1 and GX 17-1-2 (and likewise other "Z"-type low mass X-ray binaries) arise from 
jy-^ ■ the same mechanism. 

OO 
O _ 

0^ I Subject headings: accretion, accretion disks — relativity — pulsars: general — stars: neutron 

' rotation — X-rays: stars 

o 

, The variety of QPO modes from old accreting neutron stars in low mass X-ray binaries (LMXRBs) (see 

e.g. van der Klis 1995) has been recently enriched by the discovery of kiloHertz oscillations in light curves 
obtained with the Rossi X-ray Timing Explorer (RXTE). These are the first QPOs that involve timescales 
comparable to the dynamical timescales in the vicinity of a collapsed star. There are currently eleven LMXRBs 
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' from which QPOs have been detected between frequencies of ~ 0.3 and 1.2 kHz. 

Most of these are type I X-ray bursters with persistent luminosities in the 10'^^ — 10^^ ergs s~^ range, that 
according to the classification of Hasinger and van der Klis (1989) either belong to (4U 1728-34, 4U 1608-52, 
4U 1636-536, 4U 1735-444 and 4U 1820-30) or are suspected members of (4U 0614-1-091 and KS 1731-260) the 
"Atoll" group. Sco X-1, GX 5-1 and GX 17+2 are instead high luminosity "Z"-type LMXRBs. 

A pair of kHz QPO peaks is often present in the power spectra. This is the case for the persistent emission 
of the five Atoll sources, at least on occasions and/or over some range of count rates. In 4U 1728-34 and 
4U 0614-1-091 the centroid frequency of the QPOs is positively correlated with the source count rate across 
the factor of ~ 2 QPO frequency variations (from ^ 500 to ~ 1100 Hz) observed from both sources. When a 
pair of peaks is present, the difference Av between their centroid frequencies vi and 1^2 is consistent with being 
constant. In 4U 1728-34 Ai/ remains around ~ 363 Hz (within a few percent) as the frequency of the kHz 
QPOs changes by ~ 20%. A nearly coherent signal at ~ 363 Hz has been detected during a number of X-ray 
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bursts (Strohmayer et al. 1996a). In the persistent emission of 4U 0614+091 the separation Av is consistent 
with ~ 323 Hz (within ~ 10%), across a ~ 50% variation of the QPO frequency and a ~ 100 d interval between 
two RXTE observations. During part of an observation, a third narrower peak has been marginaUy detected at 
a frequency of ^ 328 Hz, consistent with the separation of the pair of kHz QPO peaks (Ford et al. 1997). In 
the persistent emission of KS 1731-260, two kHz QPO peaks have been detected around ^ 900 and ^ 1160 Hz. 
The frequency difference (~ 260 Hz) between these QPOs is consistent with half the frequency of the nearly 
coherent 524 Hz signal observed in a type I burst from this source (Smith et al. 1997; Wijnands & van der Klis 
1997). kHz QPOs have also been discovered in 4U 1608-52 (Berger et al. 1996) and 4U 1735-444 (Wijnands et 
al. 1996). 

A QPO peak with an intensity-dependent frequency from ^ 550 to ^ 800 Hz has been detected in 
4U 1820-30. Close to the upper end of this range a second QPO peak at ^ 1070 Hz is present (Smale, Zhang 
& White 1997). In the persistent emission of 4U 1636-536 a pair of QPO peaks around 840-920 Hz and 
1150-1200 Hz have been detected, the frequency separation of which is constant around Ai/ ^ 255 ± 25 Hz. 
A nearly coherent signal at ^ 581 Hz modulates the X-ray flux of this source during type I bursts (Zhang et 
aL1996, Wijnands et aZ.1997); note that half this frequency is compatible with the separation Av. Periodic 
oscillations at 589 Hz have also been observed during type I bursts from an unknown burster in the direction 
of the galactic centre, perhaps MXB 1743-29 (Strohmayer et al. 1996b). 

GX 5-1, GX 17+2 and Sco X-1, currently the only high-luminosity "Z"-type LMXRBs in the group, 
display kHz QPOs 0.3 — 1.1 kHz) with a centroid frequency which is positively correlated with the inferred 
mass accretion rate. In GX 5-1 a pair of kHz QPO peaks has been detected over part of the horizontal branch 
(HB), which mantain a constant separation of Av ~ 325 Hz. The higher frequency QPO peak of GX 5-1 
increases from V2 ~ 570 to ~ 900 Hz as the source moves from the left to the right end of the HB (van der Klis 
et al. 1996a). In GX 17+2 a kHz QPO peak is observed to drift from - 680 to 880 Hz QPO as the source 
moves down the upper part of the so-called normal branch. A higher frequency peak separated by ~ 306 Hz 
is observed close to the lower end of this range, together with an additional QPO peak around ^ 60 Hz, 
interpreted as HB oscillations (van der Klis et al. 1997a). In Sco X-1, ~ 0.5 — 1.1 kHz QPOs have been detected 
in the normal and flaring branches; their frequency correlates well both with the inferred variations of the 
mass accretion rate and the frequency variations of the 6 — 20 Hz QPO that occur along the two branches. 
A pair of kHz QPO peaks has been detected on occasions. However their separation decreased from Av ^ 310 
to ^ 230 Hz as the higher frequency peak moved from ^ 875 to ^ 1085 Hz. Correspondingly additional QPOs 
at ^ 40 — 50 Hz were present, which are intepreted as HB QPOs (van der Klis et al. 1996b, 1997b). 

Modelling of kHz QPOs has focussed on the presence of peaks at vi and V2, the constancy of their 
separation Av and the coherence of the signal (either during type I bursts or in the persistent emission) at 
frequency Av (or 2Aj/)0. 

The application of beat frequency concepts to these signals is natural. In the magnetospheric beat 
frequency model, MBFM, originally developed to interpret the 5-50 Hz QPOs from high luminosity "Z"-type 
LMXRBs (Alpar & Shaham 1985; Lamb et al. 1985), matter inhomogeneities in the accretion disk cause 
accretion to be modulated at the beat frequency between the Keplerian frequency, uk, of the innermost region 
of the accretion disk (ri„) at the magnetospheric boundary, and the neutron star spin frequency, Vg- As r,>i 
decreases with increasing accretion rates, vk increases and so does the QPO beat frequency viy = vk — Vs- 
Of the various QPO modes with frequencies < 100 Hz in "Z"-type LMXRBs, only HB QPOs are currently 
interpreted in terms of the MBFM (for a review see Lamb 1991). The basic physical picture of the MBFM 
has been confirmed mainly through the study of QPOs in a few accreting X-ray pulsars (EXO 2030+375 
and A 0535+26 in particular), where, unlike "Z"-type LMXRBs, the neutron star spin frequency Vs is (by 
definition) clearly detected, and independent estimates of the neutron star magnetic field exist (Angelini et al. 
1989; Finger et al. 1996). 



^ Though these properties are not observed in all LMXRBs with kHz QPOs, they appear to be common, with Sco X-1 the only 
clear exception to the constancy of Av (White & Zhang 1997) 



Besides the signals at Vg and an additional QPO signal at vk is expected in the MBFM from 

inhomogeneities at the inner disk boundary. For self-luminous blobs, the amplitude of this signal should be 
higher for smaller values of ri„ , due to the larger local gravitational energy and Keplerian velocities (causing 
a larger Doppler boosting around the orbit). Similarly, a small value of Tin increases the range of inclination 
angles for which a signal at vk can be caused by the neutron star occulting the blobs at Tin or conversely, by 
the blobs at rj„ obscuring part of the radiation produced closer to the neutron star surface. 

In the MBFM for kHz QPOs from LMXRBs, the lower frequency QPOs are interpreted as the beat 
frequency {ui ~ Vb) and the higher frequency QPOs as the Keplerian frequency (f2 — vk) at the magnetospheric 
boundary. The nearly coherent signal, occasionally detected at a frequency compatible with the separation 

of the QPO pair (or with 2Ai^), provides the first direct evidence of the neutron star spin frequency Vg in 
LMXRBs containing an old neutron star (Strohmayer et al. 1996a; Ford et al. 1997). Inferred magnetospheric 
radii ^ nn = {GMY/^{2t:vk)~'^'^ = 15Mq^^v]^1^^ km, (here M = MqMq is the neutron star mass and 
uk = lO^z/ifs Hz) correspond to a couple of neutron star radii at the most. The neutron star spin frequency of 
several hundred Hz and the inferred magnetic field strength of ~ 10^ — 10^ G agree with spin-up by accretion 
torques, and scenarios for the formation of millisecond pulsars. 

Miller et al. (1996) proposed an alternative beat frequency model based on the idea that the magnetic 
field of the neutron star is so low that it does not affect the dynamics of the innermost disk regions and, yet, 
high enough to give rise to beamed radiation at the neutron star surface. We do not discuss the details of 
this model, but notice that here too the QPO signal at the highest frequency 1^2 corresponds to the innermost 
Keplerian frequency of the disk, which in this model ends at the sonic radius where the radial velocity becomes 
supersonic (as opposed to the magnetospheric radius), and the lower frequency vi is the beat frequency 
between the Keplerian frequency and the neutron star spin frequency, exactly like in the MBFM. 



2. Lense-Thirring Precession in LMXRBs with kHz QPO 



The parameters inferred from beat frequency models for the kHz QPO from LMXRBs are such that frame 
dragging, one of the general relativistic effects of the "gravitomagnetic" field associated with rotating bodies, 
is substantial and can lead to observable phenomena. In particular the orbital plane of a test particle, if not 
coaligned with the equatorial plane, must precess in a prograde fashion around the angular momentum axis of 
the rotating object. In the weak field limit the nodal precession frequency, first calculated by Lense & Thirring 
(1918) (see also Wilkins 1972), is given by vlt = GMa/ {-Kc^r^), where M and a are the mass and specific 
angular momentum of the rotating object and r is the radius of the test particle (circular) orbit. By using 
Ma = 2'kIvs, and r = {GMy/^ this can be written as 

i^LT = ^^^^^^ = 13.2 J45Mo-ll/^3l/,2.5 Hz (1) 

where / = 10'*"'' 145 g cm^ is the moment of inertia of the neutron star and Vg = 3001^52.5 Hz its spin frequency. 
By using the values of and vk inferred from the application of beat frequency models to the kHz QPOs from 
LMXRBs, it is apparent that Lense-Thirring frequencies of a few tens of Hz are expected in the innermost 
disk region which produces the QPO signal at uk- Note also that in this regime, higher order terms in a/M 
would provide a correction to Eq. (1) at a level of less than one percent. 

Since sources with a pair of kHz QPO peaks allow identification of both and the only parameter 
in Eq. (1) that is not immediately identified from observations is I/M. The theoretical uncertainty in its 
determination depends on the choice of the equation of state (EOS) and the star mass: from models of rotating 
neutron stars (Friedman, Ipser and Parker 1986; Cook, Shapiro & Teukolsky 1994), it can be seen that I/M 

varies in the range 0.5 < /45/Mo < 2 for any mass and EOS. We notice that jjlt in Eq. (1) is thus independent 
of any other uncertainty related to LMXRBs, such as source distance, mass accretion rate or efficiency of 
conversion of gravitational energy into radiative energy. 



= .2,. ' (2) 



The neutron star spin frequencies discussed here lead to a stellar oblateness which, in turn, causes the 
appearance of a quadrupole term in the gravitational potential. It is well-known that the quadrupole term 
causes a precession of orbits tilted away from the star equatorial plane, with the precession frequency oc r~'^/^ 
(as opposed to for the Lense-Thirring precession). For the small radii of interest to us, it is important to 
check that this classical precession rate is still smaller than the Lense-Thirring rate, to make sure that the 
effect we predict is a genuine relativistic one. The classical precession frequency i>ci is given by (Kahn and 
Woltjer 1959) 

3 02 cos P 
Stt^ r'^VK 

where /3 is the tilt angle off the equatorial plane, and 02 is the coefficient of the quadrupole term in an 
expansion of the gravitational potential in spherical harmonics, </> w GM/r + 02^2 (cos 0), with P2 the second 
order Legendre polynomial. Outside the star, we have 02 = GA/r^, where the coefficient A is given, in terms 
of the star moment of inertia, hy A = [I^x + lyy ~ '^Izz)l'^ , where z is the star rotation axis. For rotationally 
flattened stars, far from breakup, Ixx = lyy For rotationally flattened stars, A < 0, and both Ixx and Izz are 
increased by rotation. In the existing literature we have been able to find only the numerical values of for 
non-rotating neutron star models and the corresponding rotation-increased values, (Friedman, Ipser and 
Parker 1986, Cook, Shapiro and Teukolsky 1994). Using the rotating polytropic models calculated by James 
(1964), we derived the approximate relation I^J^ — ~ 2(/^ — /^). Assuming that this holds also for 
neutron star models, we have A ~ {I^J^ — Izz)/'^- This is probably a fairly accurate approximation because the 
rotation frequencies we are interested in, < 500 Hz, correspond to still slowly rotating neutron star models: 
in fact, perusal of the tabulated models in the references given above shows that the ratio t of rotational to 
binding energy for these Vg is ~ 0.01 — 0.04, while models nearly at breakup have t ^ 0.11, for all EOS. We 
find, defining A = —rjvl2 where / — I^J^, that the parameter 77 is roughly independent of spin frequency in 
the slow-rotation regime of interest here, and it depends, albeit weakly, upon model mass and EOS (range of 
■q « (0.5 - 3.2) X 10-2). Eventually, we find 

= -4.7/45M0" ( cos /? 4/3^2^.5 Hz . (3) 

Note that the negative sign reflects the fact that classical precession is retrograde for /3 < 7r/2. Comparison 
of this with equation |^ shows that, even at these small radii, Lense-Thirring precession is faster. The total 
precession frequency Up — vlt + is thus somewhat decreased by oblateness-induced precession. 

Modulation of the X-ray flux at the precession frequency of the inner disk boundary is expected to 
occur through the same geometrical effects producing the modulation at the Keplerian frequency. In fact, 
the precessional motion of the blobs generating the signal at vk would give rise to a signal at the precession 
frequency, both for self-luminous blobs (the projected velocity along the line of sight and, therefore, the Doppler 
boosting will be modulated at the precession frequency) and occulting, or occulted, blobs (the duration of the 
occultations will depend on the phase of the precessional motion). Moreover we note that, unlike the signal 
at UK, an X-ray modulation at the precession frequency can in principle be generated, even in the absence of 
blobs, by the changing aspect of the innermost disk ring. This would modulate the ring emission as seen from 
the earth or could partially occult the emitting regions close to the neutron star surface. 

Lense-Thirring precession of viscous accretion disks is generally neglected because, following the pioneering 
work by Bardeen and Patterson (1975), many authors found that viscosity leads the disk to lie in the equatorial 
plane of the rotating object. However, this conclusion has recently been challenged by (some of) the same 
authors: even if the innermost disk regions are confined to the equatorial plane of the neutron star, Pringle 
(1996) and Maloney, Begelman and Pringle (1996) have shown that radiation pressure can lead to corrugation 
of the disk and its lifting off the neutron star equatorial plane. Pringle (1996) has shown that even initially 
planar disks are unstable to warping, so that tilting should be generic in centrally illuminated sources. 
Precessing disks have also been recently revisited by Papaloizou et aZ.(1997). It is possible that other effects, 
such as parametric amplification (Binney 1978) of the z-oscillations of diamagnetic blobs by a dipolar magnetic 
field tilted with respect to the star spin axis, may intervene (Vietri and Stella, in preparation). 



The power spectra of 4U 1728-34 given in Fig. la clearly show a broad peak centered around a frequency 



of ^ 20, 26 and 35 Hz when the higher frequency QPO peak is at a frequency of 1^2 ^ 900, 980 and 1100 Hz, 
respectively. The width of the 20 — 35 Hz peak is of 15 — 20 Hz. The power spectrum of 4U 0614+091 
(see Fig. lb) and KS 1731-260 (see Fig. Ic) also shows a similarly broad peak centered around ~ 22 Hz and 
~ 27 Hz, when 1^2 ~ 900 Hz and 1^2 ^ 1200 Hz, respectively. We propose that these tens of Hz QPO peaks 
arise from the precession of the innermost disk regions, where the signals at the Keplerian and beat frequencies 
also originate. The thick bars in Fig. 1 represent the precession frequency calculated from the sum of the 
Lense-Thirring and classical frequencies (Eqs. 1 and 3) by using = V2 (see above), and = 363, 323 and 
262 Hz in the case of 4U 1728-34, 4U 0614-^091 and KS 1731-260, respectively In all three cases, we adopted 
the same neutron star model with EOS L, Afg = 1.97, /45/A/o — 1-98 and r\ = 10^^ (Friedman, Ipser & Parker 
1986). Lense-Thirring precession dominates, with the classical term providing a ~ 20% correction. The fact 
that the (absolute) width of the tens of Hz QPO peaks is comparable to that of the kHz QPO peaks indicates 
that lifetime broadening dominates the coherence in both cases. 

Both vk and Vs are independently determined in these three sources. Thus the range of precession 
frequencies allowed by the model is determined only by the neutron star parameters I45/M0 and, to a lesser 
extent, Mq and rj. The intervals marked by the filled squares and thick bars in Fig. 1 represent the whole 
range of precession frequencies that has been obtained by using virtually all neutron star models given by 
Friedman, Ipser & Parker (1986) and Cook, Shapiro & Teukolsky (1994) (we excluded models with Mq ^ 0.5 
and "supramassive" models). We regard the very good agreement with the frequency of the corresponding 
power spectrum peaks as a strong indication in favor of our relativistic plus classical precession interpretation. 
We note that for all three sources the QPO peaks are very close to the high-frequency end of the intervals 
corresponding to the relatively stiff EOS L and Mq ~ 2. As the frequency of the kHz QPOs from 4U 1728-34 
varied, the best fit precession frequency remained in good agreement with the observed peak, showing that the 
power law dependence (Eq. |^) is probably in the correct range. 

The detection of kHz QPOs in the HB of GX 5-1 is especially important for the modelling of QPOs 
in "Z"-type LMXRBs. If one retains the standard interpretation in which the 15 — 50 Hz QPOs in the HB 
arise from the MBFM, then a magnetosphere extending up to > 5 neutron star radii is required (implying a 
magnetic field of ~ 10^ G) which argues against any interpretation of the kHz QPOs involving Keplerian orbits 
at ~ 15 — 20 km. If on the contrary, a Keplerian beat frequency model is adopted for the kHz QPOs, then 
the 15 — 50 Hz QPOs are unlikely to arise from the MBFM. We regard the presence of a pair of kHz QPOs 
and the constancy of their frequency difference as strong evidence in favor of the latter interpretation. In this 
context, it appears natural to interpret also the HB QPOs of GX 5-1 in terms of our precession model. For 
Ps = 325 Hz and vk = 1^2 ^ 570 — 900 Hz, the maximum precession frequency (again the one from EOS 

L and Mq ~ 2) is ^ 8 — 20 Hz. Similarly in the case of GX 17+2, we obtain a maximum precession frequency 
of ~ 22 Hz, when the source is in the upper normal branch with vk — V2 988 Hz and — lS.v ^ 306 Hz, 
and HB QPOs at ~ 60 Hz are simultaneously present. These precession frequencies are clearly a factor of 
2-2.5 smaller than those of the HB QPOs. One possibility is that the neutron stars of "Z"-type LMXRBs 
have /45/Afo 4 — 5, i.e. a somewhat higher value than Friedman, Ipser & Parker (1986) and Cook, Shapiro 
& Teukolsky (1994) calculate. Alternatively the HB QPOs might arise from the second harmonics of the 
precession frequency, due to some yet unknown reason. We note that in the case of GX 5-1 the dependence of 
the HB QPO frequency on V2 is in the same range as the dependence of on vk in our model. 

In summary, we propose that the peaks observed in the power spectra of several "Atoll" LMXRBs at a 
frequency of tens of Hz (and perhaps also the HB QPOs of "Z"-type LMXRBs) are due to precession of the 
innermost disk region, which is dominated by the Lense-Thirring effect. The precession frequency depends 
upon the Keplerian frequency at the innermost disk radius, and the neutron star spin frequency, both of which 
are independently determined through the application of beat frequency models to the kHz QPOs from these 
sources. The best fit is obtained for a stiff EOS and A/q ~ 2, a value supported also from the interpretation of 
the highest observed values of V2 in terms of the Keplerian frequency of the marginally stable orbit (Kaaret, 
Ford & Chen 1997; Zhang, Strohmayer & Swank 1997). 

LS acknowledges useful discussions with D. Bini, I. Ciufolini, D. Psaltis. This work was partially supported 
through ASI grants. 
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Fig. 1.— Power spectra of the RXTE light curves of 4U 1728-34 (panel a), 4U 0614+091 (panel b) and KS 1731- 
260 (panel c) (after Strohmayer et al. 1996a, Ford et al. 1997 and Wijnands & van der Khs 1997). The thick 

bars indicate the best precession frequency obtained from the model discussed in the text. The range of allowed 
precession frequencies for a variety of neutron star EOS and masses is marked by thick bars and filled squares. 



